How to determine Fermi vectors by angle resolved photoemission? 
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Angle resolved photoemission spectroscopy (ARPES) has been commonly applied to evaluate the 
shape of Fermi surfaces by employing simple criteria for the determination of the Fermi vector kF 
parallel to the surface such as maximum photoemission intensity at the Fermi level or discontinuity 
in the momentum distribution function. Here we show that these criteria may lead to large uncer- 
tainties in particular for narrow band systems. We develop a reliable method for the determination 
of Fermi vectors employing high resolution ARPES at different temperatures. The relevance and 
accuracy of the method is discussed on data of the quasi two dimensional system TiTe2- 



A wide variety of physical phenomena of crystalline 
materials, such as e. g. transport, optical and magnetic 
response, and phase transitions rely on details of the 
topology of the Fermi surface (FS) . Its experimental de- 
termination performed by traditional techniques like de 
Haas- van Alphen effect, magnetoacoustic effect, Comp- 
ton scattering, or positron annihilation have, however, 
been restricted to bulk materials. All techniques in com- 
mon is the indirect information on the shape of the Fermi 
surface. While the first two determine extremal cross- 
sections of FS's in a plane normal to the applied mag- 
netic field, the latter yield information on one and two 
dimensional projections of FS's. More complex cases, 
such as superlattices, heterostructures or even clean sur- 
faces can also hardly be accessed. Angle-resolved pho- 
toemission spectroscopy (ARPES) has emerged as prob- 
ably the most powerful tool for determining the occupied 
electronic band structure of solids and their surfaces. Re- 
cently, it has been extensively applied to gain insight into 
the topology of Fermi surfaces of a variety of materials 
ranging from conventional three-dimensional metals like 
W and Cu [pl|]|] to quasi two dimensional layered com- 
pounds ^,|]6 , purple bronzes Q and high Tc cuprate 
materials [gflE]. The accuracy of the determination of 
the Fermi surface by ARPES, however, has never been 
questioned and it turns out that even for the extensively 
studied BiSrCaCuO it is not at all clear that the topology 
of the normal state FS shows hole- like pockets around the 
corners of the Brillouin zone p|,ppC| ] or electron-pockets 
around the center of the BZ 

It is widely assumed that ARPES measures the spec- 
tral function A(k, iS) of the one-particle system times the 
Fermi function f(u>) and matrix elements do not play 
a significant role (see e.g. ref. p2| , p^| ). This motivates 
studies of band dispersions, line shapes, momentum dis- 
tribution functions, and Fermi surfaces. Fermi vectors 
have been extracted from ARPES data employing cri- 
teria like (i) maximum ARPES intensity at the Fermi 
level Ep 0,@,0,||Jl4| , (ii) ma xlVfcj o f the energy integrated 
photoemission intensity [^p| ^5| , |l6|| or (iii) fitting ARPES 
peak positions over several emission angles and extrapo- 



lating the dispersion to Ep p7| . However, none of these 
techniques explicitly considers the detailed mechanism of 
the photoemission process. In particular, matrix clement 
effects and differences between photocurrcnt and spec- 
tral function have been totally neglected and make these 
simple interpretations highly questionable. Although, in 
most cases, photoemission peak positions still resemble 
locations of electronic bands, the shape of spectral func- 
tions is in general not directly reproduced by the pho- 
tocurrent []l8| , [l9|| . If photoemission calculations within 
the one-step model are not available reliable simple pro- 
cedures for an analysis of high resolution photoemission 
data are needed. 

Employing high resolution photoemission spectroscopy 
we show in this Letter how Fermi vectors can be deter- 
mined with high accuracy when intensity modifications 
due to the photoemission process are explicitly elimi- 
nated by comparing photoemission spectra taken at dif- 
ferent temperatures. The reliability and applicability of 
the method will be demonstrated on the layered Fermi 
liquid reference material IT — TiTe2. It does not show 
an indication of charge density waves and exhibits Ti 3d 
conduction band emissions which are well separated from 
other bands @|I]j22). It may thus serve as a model sys- 
tem for assessing the accuracy of the determination of 
the Fermi surface of quasi two dimensional systems. 

All results described here were obtained on clean chem- 
ical vapor transport grown IT — TiTe2 samples prepared 
by cleavage in ultra high vacuum. Photoemission spectra 
were taken with synchrotron radiation supplied from the 
storage ring DORIS III at Hamburg Synchrotron Radi- 
ation Laboratory (HASYLAB) using our Angular Spec- 
trometer for Photoelectrons with High Energy REsolu- 
tion (ASPHERE). The electrons were detected by use of a 
180° spherical analyzer mounted on a two-axes goniome- 
ter with an absolute angular precision of better than 0.1°. 
The overall energy resolution was chosen to 30 meV. The 
position of the Fermi level was determined from spectra 
of polycrystalline gold with an accuracy of ±1 meV. 

The photocurrent within the sudden approximation is 
given by the spectral representation of the one-particle 
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Green function: Z(k, lo) = Io(k)f(Lu)A(k, lo) where f is 
the Fermi function and A(k, lu) is the one-particle spec- 
tral function. The prefactor Iq involves the transition 
matrix element and is thus k dependent. As will be 
shown below this k dependence crucially affects the reli- 
ability of the determination of kp . Depending on photon 
energy the component of k perpendicular to the surface 
is varied in angle resolved photoemission. It further 
complicates a full three dimensional determination of the 
Fermi vector. This, however, is going beyond the scope 
of this Letter and will be discussed in a forthcoming pub- 
lication. 

The component of the Fermi vector parallel to the sur- 
face kp|| can be identified by locating the peak of A(k, lu) 
at zero binding energy. This is not a trivial matter 
since the photocurrent corresponds to /(w)A(k, lu) and 
the peak in the photocurrent is what results when the 
Fermi function cuts off the peak of the spectral function 
A(k, lo). Different approaches have been used to solve 
this difficulty: (i) find the peak in 7(k, Ep) (maximum in- 
tensity at the Fermi level) , (ii) use the sum rule [^3| relat- 
ing the energy integrated spectral function to the momen- 
tum distribution function: n(k) = J_ dto f (lu) A(k., lu). 
By definition the momentum vectors kp constituting the 
Fermi surface are identified through a step or at least 
the vertical slope of the momentum electron distribution 
n(k) at zero temperature [p|Jl5|. Method (iii) employs a 
fit of photoemission peak positions over a wide range of 
emission angles and extrapolation to zero binding energy 
to obtain an estimate for kp||. 

In Fig. (a) and (b) we show the results of a simulation 
for a narrow band represented by a gaussian peak (width: 
AE F whm = 100 meV, dispersion: AE/Ak = -0.2 eVA) 
crossing the Fermi level at k|| = OA -1 . Fermi vectors de- 
termined using the approaches (i) and (ii) are plotted as a 
function of the experimental resolution. Solid lines show 
results for a matrix element Io(k) — const, shaded areas 
for a weakly changing Io(k) ~ a(k + 1), — 1 < a < +1. 
Both approaches provide values for k F || revealing system- 
atic deviations of up to about 20% (with respect to the 
distance TM = 0.96A -1 of the Brillouin zone) strongly 
depending on the experimental resolution. It should be 
noted here that these approaches may be additionally af- 
fected in strongly correlated Fermi liquid systems where 
incoherent backgrounds (e.g. electron hole pairs) domi- 
nate the modulation of n(k) |24 ], 

A reliable approach for the determination of kp|| is to 
compare spectra at different temperatures (AT method). 
Centered at Ep and kp and with respect to u>, the differ- 
ence of the Fermi functions for temperatures Ti and T2 
is an odd function, A(k.p,Lu) is even, and therefore the 
integral of their product over a symmetric energy win- 
dow vanishes. In practice, the integration is replaced by 
the finite energy resolution given by a symmetric analyzer 
function w, which vanishes outside [— e, e]. The difference 
of the intensities taken at Ep = is 



A/(k||) = / (k||) J <b>AQL,u)[f Tl (L>) ~ /t»Mw) , 

assuming that for the intervals of interest A is indepen- 
dent of temperature and Iq of temperature and energy. 
It then follows from A/(k||) = that k|| = kpi|. This 
is valid for all A which satisfy in [— e, e] that A(k.p,io) — 
A(k F , -lu) and A(k, u) ^ A(k, -lo) for k ^ k F . Then at 
general k the differences of A do not change sign and 

A/ = Jo J dLu(A(k,Lu) - A(k, -Lu))[f Tl - /t> ? . 

A variety of spectral functions fulfill these conditions, in- 
cluding Lorentzians, Gaussians, Voigt profiles, the Lut- 
tinger model iS, an extension by Matho |2l]j2(J, two 
dimensional [2j and marginal pi| Fermi liquids. A high 
energy resolution extends the scope. The independence 
of temperature of such integrals in the vicinity of kp and 
Ep has already been tested [H) and will later be shown 
again for the present case. With a high k|| resolution 
the criterion is still valid since for a symmetric band dis- 
persion (with respect to (Ep, kp)) the additional contri- 
butions annihilate. The accuracy might be lowered due 
to k dependencies of Iq. This possible influence, how- 
ever, can directly be excluded when the experimental k\\ 
resolution windows lies well within a regime where point 
symmetry of A I(k\\ ) around AI = is observed in the ex- 
perimental data. Note that for a tiny energy interval the 
required fixed-k mode of photoemission coincides with 
the employed fixed-angle mode well within the angular 
resolution. A simulation of the AT method for a gaus- 
sian profile and the Luttinger model is shown in Fig. [I] (c) 
and (d). The Fermi vector k|| — is exactly determined 
by AI = in both cases for all resolutions. 

In Fig. ^| (a) we show high energy and angle resolved 
photoemission spectra of IT — TiTe2 along the TM di- 
rection of the Brillouin zone. For the Ti 3d-band which 
is well separated from contributions of other bands one 
observes besides the crossing of the band connected with 
an intensity break down a characteristic behavior of the 
line shape. Line shape studies have been performed em- 
ploying either a Fermi liquid scenario using a many body 
spectral function with an imaginary part of the self- 
energy depending quadratically on energy referred to Ep 
pof or profiles suggested by Matho |2l|]26H for approxi- 
mating the spectral function A(k, lo). In the close vicin- 
ity of the Fermi vector the experimental spectra could 
excellently be fitted. But for emission angles larger than 
~ 18° the quality of the fits is drastically reduced. Values 
obtained for the Fermi vector were varying between 14° 
(0.51A- 1 ) §l| and 14.75° (0.53A" 1 ) @||,||. These 
studies, however, did not consider the photoemission pro- 
cess explicitly and intrinsic spectral information can be 
distorted by varying matrix elements or other secondary 
effects. This problem is further illustrated in Fig. |^(a) 
where an excellent fit (solid lines) can already be achieved 
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for all emission angles employing simple Voigt profiles 
times the Fermi function convoluted by the spectrome- 
ter response. This demonstrates that without explicit 
consideration of the photoemission process sophisticated 
spectral functions can hardly be delineated. Peak posi- 
tions resulting from the fit give the experimental disper- 
sion shown in Fig. |2| (b). The band crossing is observed at 
16.7° corresponding to 0.60 A" 1 . It should be noted here 
that in addition to an increased experimental and calcu- 
lational (fit) effort the values for kpii obtained by method 
(iii) rely on the profiles employed. For the Matho profile 
the band crossing is observed at 14.75° (see Fig. || (b)). 

For an application of the AT method we plot normal- 
ized photoemission spectra of IT — TiTe2 taken at 100 K 
and 190 K (Fig. || (a)). The broadening of the Fermi 
function according to temperature is evident. Intensities 
at the Fermi level (AE = 30 meV) and intensities inte- 
grated over the whole spectrum are depicted in (b) and 
(c) , respectively. According to the AT method the Fermi 
vector is given by the intersection of the curves in (b) 
which can clearly be identified at a value of 16.6° ± 0.1° 
(0.598±0.004A^ 1 ). For comparison we also show the val- 
ues for kp obtained using the maximum intensity (i) and 
maximum gradient (ii) methods described above. The 
maximum intensity at Ep can be observed on a relatively 
broad peak at around 17.8° (0.64A -1 ) showing the sys- 
tematic erroneous shift of k F || towards occupied states as 
already observed in the simulation of Fig. |l]. The energy 
integrated intensity and its derivative are depicted in Fig. 
H (c) . Since the gradient of the integrated intensity is only 
marginally changing in the k|| regime of interest a rather 
broad maximum is observed making a detailed quantita- 
tive comparison with FS calculations very difficult. Sym- 
metrized spectra /(kp, u>) + /(kp, — uA = Io(kF)A(k-F,ui) 
eliminating the Fermi function |29l , |l^| are shown in Fig. 
| (d). For the spectra taken at 100 K and 190 K no dif- 
ferences can be observed ruling out any T-dependence of 
io(kF)^4(kF, uj) in this temperature range. 

In summary, we have developed a reliable and sim- 
ple method to determine Fermi vectors by angle resolved 
photoemission spectroscopy. Employing temperature dif- 
ference spectra the photoemission process hiding the 
spectral function in measured spectra is explicitly con- 
sidered. Therefore, the method even works for systems 
showing rapidly changing matrix elements with variation 
of k. For the layered material IT — TiTe2 we have demon- 
strated that the accuracy of the determined kp values can 
be significantly improved to better than ±0.4% of the di- 
mension of the Brillouin zone. Compared to other meth- 
ods revealing up to one order of magnitude larger error 
bars together with systematic deviations of up to 20% re- 
liable quantitative comparisons between experiment and 
Fermi surface calculations going beyond similarities of 
the shapes will now become possible. 

Experimental help by R. Schwedhelm and S. Woedtke 
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FIG. 1. Systematic deviations for k F || determined by ex- 
perimental procedures employing methods (i) max(J(k, Ef)) 
and (ii) max(|Vfc7i(fc)|) plotted as a function of experimental 
resolution (gaussian type) for 300 K (a) and 30 K (b). Sim- 
ulation for a gaussian peak (width: AEfwhm = 100 meV, 
dispersion: AE/Ak = -0.2 eVA) crossing the Fermi level at 
fc F || = A" 1 . Solid lines and shaded areas represent simula- 
tions for constant and varying matrix elements, respectively. 
Simulation results for the AT method are shown in (c) and 
(d). See text for details. 

FIG. 2. (a) Angle resolved photoemission spectra (EDC) 
associated with the Ti 3d-band of IT — TiTe2 in the direction 
TM taken with high resolution and line shape fit using Voigt 
profiles, (b) Band dispersion obtained from the fitted peak 
positions (method (hi)) using Voigt and Matho profiles. Note 
the different band crossings at Ef. 

FIG. 3. (a) Angle resolved photoemission spectra (AE = 
30 meV) associated with the Ti 3d-band of IT - TiTe 2 in 
the direction TM taken at 100 K and 190 K. (b) Intensi- 
ties at the Fermi level (method (i)). The intersection marks 
the Fermi vector (AT method). An uncertainty of ±0.1° 
(hatched area) emerges from the error bar of the Fermi energy 
(±1 meV) . (c) Intensities integrated over the whole spec- 
trum and k-derivative (smoothed) of the integrated intensity 
(method (ii)). Grey bars mark the uncertainties for the de- 
termination of k F (95% of the maximum value), (d) Sym- 
metrized spectra showing no differences for spectral functions 
at 100 K and 190 K. 
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